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Soufriere  Hills  Volcano  had  two  periods  of  repetitive  Vulcanian  activity  in  1997.  Each  explosion  discharged 
the  contents  of  the  upper  0.5-2  km  of  the  conduit  as  pyroclastic  flows  and  fallout:  frothy  pumices  from  a 
deep,  gas-rich  zone,  lava  and  breadcrust  bombs  from  a  degassed  lava  plug,  and  dense  pumices  from  a 
transition  zone.  Vesicles  constitute  1-66  vol.%  of  breadcrust  bombs  and  24-79%  of  pumices,  all  those  larger 
than  a  few  tens  of  pm  being  interconnected.  Small  vesicles  (<few  tens  of  pm)  in  all  pyroclasts  are  interpreted 
as  having  formed  syn-explosively,  as  shown  by  their  presence  in  breadcrust  bombs  formed  from  originally 
non-vesicular  magma.  Most  large  vesicles  (>few  hundreds  of  pm)  in  pumices  are  interpreted  as  pre-dating 
explosion,  implying  pre-explosive  conduit  porosities  up  to  55%.  About  a  sixth  of  large  vesicles  in  pumices, 
and  all  those  in  breadcrust  bombs,  are  angular  voids  formed  by  syn-explosive  fracturing  of  amphibole 
phenocrysts.  An  intermediate-sized  vesicle  population  formed  by  coalescence  of  the  small  syn-explosive 
bubbles.  Bubble  nucleation  took  place  heterogeneously  on  titanomagnetite,  number  densities  of  which 
greatly  exceed  those  of  vesicles,  and  growth  took  place  mainly  by  decompression. 

Development  of  pyroclast  vesicle  textures  was  controlled  by  the  time  interval  between  the  onset  of  explosion- 
decompression  and  surface  quench  in  contact  with  air.  Lava-plug  fragments  entered  the  air  quickly  after 
fragmentation  (~  1 0  s),  so  the  interiors  continued  to  vesiculate  once  the  rinds  had  quenched,  forming  breadcrust 
bombs.  Deeper,  gas-rich  magma  took  longer  (~50  s)  to  reach  the  surface,  and  vesiculation  of  resulting  pumice 
clasts  was  essentially  complete  prior  to  surface  quench.  This  accounts  for  the  absence  of  breadcrusting  on  pumice 
clasts,  and  for  the  textural  similarity  between  pyroclastic  flow  and  fallout  pumices,  despite  different  thermal 
histories  after  leaving  the  vent.  It  also  allowed  syn-explosive  coalescence  to  proceed  further  in  the  pumices  than 
in  the  breadcrust  bombs.  Uniaxial  boudinage  of  amphibole  phenocrysts  in  pumices  implies  significant  syn- 
explosive  vesiculation  even  prior  to  magma  fragmentation,  probably  in  a  zone  of  steep  pressure  gradient  beneath 
the  descending  fragmentation  front.  Syn-explosive  decompression  rates  estimated  from  vesicle  number  densities 
(>0.3-6.5  MPa  s-1)  are  consistent  with  those  predicted  by  previously  published  numerical  models. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Explosive  volcanic  eruptions  are  driven  by  the  nucleation,  growth 
and  coalescence  of  gas  bubbles,  followed  by  fragmentation  of  the 
magmatic  foam  into  a  suspension  of  pyroclasts  and  gas  that  is 
discharged  at  high  velocities  into  the  atmosphere.  Studies  of  pyroclast 
textures,  coupled  with  experimental  and  numerical  approaches,  have 
advanced  understanding  of  these  processes  (Lensky  et  al.  2004; 
Spieler  et  al.  2004b;  Adams  et  al.  2006;  Toramaru  2006;  Gardner 
2007;  Cluzel  et  al.  2008;  Koyaguchi  et  al.  2008  and  references  therein), 
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but  many  questions  remain.  One  concerns  the  relative  importance  of 
homogeneous  versus  heterogeneous  nucleation.  Homogeneous  nu¬ 
cleation  requires  gas  supersaturations  of  at  least  several  tens  of  MPa 
(Mangan  and  Sisson  2000;  Mourtada-Bonnefoi  and  Laporte  2002; 
2004;  Mangan  et  al.  2004),  whereas  heterogeneous  nucleation 
requires  lower  supersaturations  (Hurwitz  and  Navon  1994;  Gardner 
2007;  Cluzel  et  al.,  2008).  The  degree  of  equilibrium  between  gas  and 
melt  during  bubble  growth  also  has  an  effect.  Equilibrium  degassing 
requires  efficient  volatile  diffusion  coupled  with  melt  viscosity  low 
enough  to  allow  free  gas  expansion  (Lyakhovsky  et  al.,  1996;  Liu  and 
Zhang,  2000;  Lensky  et  al.,  2004).  High  degrees  of  disequilibrium 
favour  short-lived  eruptions,  whereas  equilibrium  allows  more 
sustained  fragmentation  (Melnik  and  Sparks,  2002;  Mason  et  al, 
2006).  Another  issue  concerns  the  timing  of  bubble  growth  and 
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coalescence  relative  to  fragmentation  and  eruption.  Some  authors 
postulate  little  growth  following  fragmentation  (Klug  and  Cashman, 
1991)  whereas  others  envisage  significant  post-fragmentation 
growth  (Thomas  et  al.  1994;  Kaminski  and  Jaupart,  1997).  Post¬ 
fragmentation  bubble  growth  is  largely  controlled  by  melt  viscosity, 
being  important  in  mafic  melts  and  less  so  in  silicic  melts  with 
viscosities  >108-109  Pa  s  (Thomas  et  al.,  1994;  Gardner  et  al.,  1996; 
Kaminski  and  Jaupart,  1997).  Bubble  coalescence  and  connections 
control  permeability  acquisition  and  the  ability  of  magma  to  outgas 
during  ascent. 

Vesicle  size  distributions  provide  information  on  magma  vesicu- 
lation  history.  Pumices  commonly  contain  multiple  vesicle  popula¬ 
tions  covering  a  large  range  of  sizes  (Klug  and  Cashman,  1996;  Klug 
et  al.,  2002;  Adams  et  al.,  2006)  that  may  result  from  coalescence 
following  a  single  nucleation  event  (Orsi  et  al.,  1992;  Klug  and 
Cashman,  1994, 1996;  Klug  et  al.,  2002;  Burgisser  and  Gardner,  2005). 
Alternatively,  each  population  may  represent  a  distinct  nucleation 
event,  consistent  with  some  ascent  models  which  predict  multiple 
events  for  viscous  magma  (Witham  and  Sparks,  1986;  Proussevitch 
and  Sahagian,  1996;  Blower  et  al.,  2001 ;  Massol  and  Koyaguchi,  2005). 
Small  vesicles  are  commonly  attributed  to  syn-explosive  vesiculation 
that  generates  an  exponential  size  distribution  (Mangan  et  al.,  1993; 
Klug  and  Cashman,  1996;  Klug  et  al.,  2002;  Adams  et  al.,  2006).  Size 
distributions  of  larger  populations  typically  obey  power  laws  usually 
attributed  to  coalescence  (Klug  et  al.,  2002;  Houghton  et  al.,  2003; 
Gurioli  et  al.,  2005;  Adams  et  al.,  2006;  Klug  and  Cashman,  1996), 
although  multiple  nucleation  events  also  generate  power-law 
distributions  (Blower  et  al.,  2001).  Magma  decompression  rates  can 
be  estimated  from  vesicle  number  densities  assuming  a  unique  and 
brief  nucleation  event  (Toramaru,  2006;  Cluzel  et  al.,  2008). 

Detailed  studies  of  eruptive  products  are  required  to  address 
these  questions  and  provide  ground  truth  for  models.  Most  vesiculation 
studies  to  date  have  concerned  Plinian  eruptions.  In  this  paper  we  study 
vesiculation  during  a  sequence  of  well  documented  Vulcanian  explo¬ 
sions  at  Soufriere  Hills  Volcano  in  1997.  The  explosions  have  been 
previously  described  (Druitt  et  al.,  2002;  Cole  et  al.,  2002)  and  modelled 
(Melnik  and  Sparks,  2002,  Clarke  et  al.,  2002;  Formenti  et  al.,  2003;  Diller 
et  al.,  2006;  Mason  et  al.,  2006),  and  their  products  studied  texturally 
(Formenti  and  Druitt,  2003;  Clarke  et  al.,  2007)  and  chemically  (Harford 
et  al.,  2003).  A  key  feature  was  the  eruption  of  pyroclasts  of  a  wide  range 
of  types,  including  dense  lava  fragments,  breadcrust  bombs  and  pumices 
of  different  densities.  Textural  analysis,  including  a  set  of  high- 
resolution  vesicle-size  distributions,  enables  us  to  recognize  populations 
of  vesicles  formed  by  explosion  decompression,  quantify  bubble 
nucleation  mechanisms  and  decompression  rates,  and  constrain  the 
timing  of  bubble  nucleation,  growth  and  coalescence  during,  and 
immediately  following,  a  typical  explosion.  In  a  companion  paper  we 
present  measurements  of  groundmass  water  contents  and  reconstruct 
the  state  of  the  pre-explosion  conduit  (Burgisser  et  al.,  in  press). 

2.  The  1997  Vulcanian  explosions  of  the  Soufriere  Hills  Volcano 

The  eruption  of  Soufriere  Hills  Volcano  (Fig.  1 )  began  phreatically  in 
July  1995;  extrusion  of  lava  began  in  November  of  the  same  year  and 
continued  intermittently  until  the  time  of  writing.  The  explosions  in 
1997  occurred  every  3-63  h  (mean  of  ~10  h)  in  two  periods:  thirteen 
between  4  and  12  August,  and  seventy  five  between  22  September  and 
21  October  (Druitt  et  al.,  2002).  Each  consisted  of  an  initial  high- 
intensity  phase  lasting  a  few  tens  of  seconds,  followed  by  a  waning 
phase  lasting  1-3  h.  Multiple  jets  were  ejected  at  40-140  m  s~ 1  during 
the  first  10-20  s  of  each  explosion,  then  collapsed  back  to  form 
pumiceous  pyroclastic  flows  that  travelled  up  to  6  km  from  the  crater 
(Formenti  et  al.,  2003).  Fallout  of  pumice  and  ash  occurred  from  high 
(3-15  km)  buoyant  plumes  that  developed  above  the  collapsing 
fountains.  Fallout  and  flow  took  place  at  the  same  time  from  individual 
explosions.  Each  explosion  discharged  on  average  8  x  108  kg  of  magma, 


about  two-thirds  as  pyroclastic  flows  and  one-third  as  fallout, 
representing  a  conduit  drawdown  of  0.5-2  km  (Druitt  et  al.,  2002). 
Studies  of  quench  pressures  using  microlite  contents  and  glass  water 
contents  support  a  maximum  drawdown  of  ~2  km  (Clarke  et  al.,  2007; 
Burgisser  et  al.,  in  press).  Each  explosion  started  when  magma 
overpressure  exceeded  the  strength  of  an  overlying  degassed  plug  and  a 
fragmentation  front  propagated  down  the  conduit  at  a  few  tens  of  m  s_  1 
(Druitt  et  al.,  2002;  Clarke  et  al.,  2002;  Melnik  and  Sparks,  2002;  Spieler 
et  al.,  2004a;  Diller  et  al.,  2006;  Mason  et  al.,  2006).  After  each  explosion, 
magma  rose  up  the  conduit  before  the  onset  of  a  new  explosion. 

The  Soufriere  Hills  andesite  contains  phenocrysts  of  plagioclase, 
hornblende,  orthopyroxene,  magnetite,  ilmenite  and  quartz  set  in 
rhyolitic  glass.  The  pre-eruptive  temperature  was  ~  850  °C  (Devine  et  al., 
2003). 

3.  Methodolgy 

Field  work  was  carried  out  in  2006  and  2008  at  three  sites  (Fig.  1 ): 
sites  1  and  2  are  situated  on  the  fans  of  overlapping  pyroclastic  flow 
lobes  from  the  explosions,  and  site  3  is  a  composite  layer  of  fallout 
pumice  from  many  explosions.  Fallout  pumices  were  also  collected  at 
a  fourth  site  (site  4;  Fig.  1)  during  an  explosion  in  August  1997.  Field 
descriptions  were  made  using  a  rock  saw  to  cut  perpendicular  to  any 
flow  banding  and  parallel  to  any  crystal  fabric,  and  over  100 
representative  pyroclasts  were  taken  for  laboratory  study. 

Abundances  of  isolated  and  connected  vesicles  were  measured  on 
2-5  cm  cubes  cut  from  30  breadcrust  bombs  and  34  flow  and  fall 
pumices  using  a  Multivolume  1305  Helium  Pycnometer  and  the 
method  of  Formenti  and  Druitt  (2003),  which  is  explained  in  the 
Supplementary  electronic  material. 

Separate  measurements  were  made  on  the  rims  and  cores  of  22 
breadcrust  bombs.  Twenty-six  of  the  pumice  clasts  ranged  from  lapilli  to 
block  size,  all  being  <20  cm  in  diameter.  Measurements  were  also  made 
on  multiple  core-to-rim  samples  from  eight  pumices  >30  cm  in  diameter. 
Texturally  or  compositionally  banded  pyroclasts  were  not  included. 

Microscopic  observations  were  made  on  the  broken  surfaces  of 
pyroclast  fragments  using  a  Jeol  JSM-591LV  Scanning  Electron  Micro¬ 
scope  (SEM)  at  an  acceleration  of  15kV,  and  on  polished  epoxy- 
impregnated  thin  sections  using  the  SEM  and  a  stereomicroscope. 

Six  samples  representative  (in  terms  of  vesicularity  and  texture)  of 
the  pyroclast  assemblage  were  chosen  for  high-resolution  analysis 
of  vesicle  and  crystal  size  distributions.  Banded  clasts,  and  those  with  a 
significant  fraction  of  non-spherical  vesicles,  were  excluded,  thereby 
justifying  use  of  a  single,  randomly  oriented  thin  section  for  each 
sample.  Vesicle  and  crystal  size  distributions  were  measured  by  image 
analysis  in  two  dimensions  (Toramaru,  1990;  Mangan  et  al.,  1993; 
Klug  and  Cashman  1 994, 1 996;  Klug  et  al.,  2002;  Adams  et  al.,  2006;  Shea 
et  al.,  2010).  The  technique,  described  fully  in  Appendix  A,  allowed 
objects  as  small  as  ~  1  pm  to  be  measured.  Differential  epoxy  penetration 
enabled  us  to  distinguish  interconnected  from  isolated  vesicles.  To 
represent  the  state  of  the  magma  immediately  prior  to  the  last  discernible 
stage  of  coalescence,  we  manually  ‘decoalesced’  neighbouring  vesicles 
separated  by  a  partially  retracted  wall.  Volume  distributions  were 
assumed  to  equal  area  distributions  (Klug  et  al.,  2002).  Volumetric 
number  densities  (Nv)  were  calculated  from  area  number  densities  (Na) 
using  both  the  methods  of  Cheng  and  Lemlich  ( 1 983 )  and  of  Sahagian  and 
Proussevitch  (1998),  which  yield  very  similar  values  (Table  1 ).  Values  of 
Nv  presented  in  this  paper  are  those  obtained  using  the  first  method,  for 
reasons  discussed  in  the  Appendix  A. 

4.  Field  descriptions 

The  pyroclast  assemblage  consists  predominantly  of  pumices  of 
different  colours,  vesicularities  and  textures,  with  less  than  a  few 
percent  of  breadcrust  bombs  and  dense  glassy  lava  clasts.  Pyroclasts  of 
all  types  were  present  in  the  pyroclastic  flow  deposits,  although  the 
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Fig.  1.  Map  of  Montserrat  showing  Soufriere  Hills  Volcano  and  sampling  locations.  Grey:  pyroclastic  flow  deposits  of  the  1997  Vulcanian  explosions. 


relative  proportions  varied  from  lobe  to  lobe,  while  dense  lava  and 
breadcrust  bombs  were  absent  in  the  fallout.  The  samples  described 
below  come  from  several  different  explosions,  and  cannot  be  assigned 
to  specific  dates/times  owing  to  the  complex  superposition  of  flow 
and  fallout  lobes  from  the  many  events.  They  represent  the  products 
of  an  ‘average’  explosion,  as  justified  by  (1 )  the  first-order  similarity  of 
all  the  explosions  (Druitt  et  al.,  2002),  and  (2)  the  presence  of  the 
entire  pyroclast  spectrum  in  all  pyroclastic  flow  lobes  we  examined. 

Pumices  in  the  pyroclastic-flow  deposits  occur  as  lapilli  and  blocks 
up  to  >1  m  in  diameter  with  subangular-to-rounded  shapes  due  to 
abrasion  during  transport.  They  range  from  beige,  well  vesiculated 
varieties,  to  grey,  brown  or  black  denser  varieties  (Fig.  2a-b).  A  pink 
colouration  affects  the  surfaces  of  many  blocks,  but  rarely  pervades 
the  interiors.  While  the  majority  of  pumices  are  texturally  homoge¬ 
neous  in  hand  specimen,  some  denser  ones  are  flow  banded  with 
phenocryst  alignment  in  the  plane  of  banding.  Rare  compositional 
banding  defined  by  trails  of  disintegrated  mafic  inclusions  also  occurs. 
All  pumice  clasts  (as  distinguished  from  breadcrust  bombs)  lack 
surface  breadcrusting.  This  probably  cannot  be  explained  by  abrasion, 
because  breadcrust  fragments  are  not  observed  in  the  flow  matrices. 
All  pumices  smaller  than  ~30cm  lack  radial  gradients  in  vesicle 
abundance  or  size.  However,  some  blocks  larger  than  this  exhibit 
visibly  obvious  radial  gradients  in  vesicle  size,  with  an  outer  3-7-cm- 
thick  rind  with  vesicles  up  to  several  mm,  and  a  more  coarsely 
vesicular  interior  containing  vesicles  up  to  an  order  of  magnitude 
larger  (Fig.  2c).  In  some  cases  a  crude  cm-scale  radial  jointing  affects 
the  rind.  The  rind  is  inferred  to  represent  the  initial  textural  state  of 
the  pumice,  while  the  interior  records  vesicle  coarsening  that  took 
place  during  or  after  emplacement.  The  possibility  that  the  interior 
represents  the  initial  state,  and  that  the  rind  developed  by  compaction 
during  rolling  in  the  pyroclastic  flow,  is  not  favoured  because  (1)  the 
rinds  texturally  resemble  the  majority  of  smaller  pumice  blocks  and 
lapilli,  whereas  vesicles  in  the  interiors  are  abnormally  coarse,  and 
(2)  no  circumferential  flattening  of  rind  vesicles  is  observed.  Many 
blocks  also  contain  large  voids  up  to  several  cm  across,  including 


anastamosing  vesicle  pipes  and  channels,  ductile  tears  in  the  plane  of 
flow  banding,  and  curviplanar  tears  and  cracks  subparallel  to  clast 
margins  (Fig.  2d),  which  together  account  for  <10%  of  the  total 
vesicularity. 

Fallout  pumices  are  up  to  several  cm  in  size  and  most  preserve 
their  original  eruption-fragmentation  shapes,  unmodified  by  abrasion 
in  pyroclastic  flows  or  breakage  on  ground  impact.  They  range  in 
colour  from  white  to  brown  and  in  shape  from  spheroidal  to  tabular, 
the  latter  comprising  about  three-quarters  of  the  sample  suite.  Again, 
no  surface  breadcrusting  is  observed. 

Breadcrust  bombs  occur  from  a  few  cm  to  over  a  metre  in  diameter. 
They  have  vesicular  interiors  surrounded  by  darker,  less  vesicular 
<10  mm  glassy  rinds.  A  continuous  range  of  textural  varieties  are 
observed  between  two  endmembers.  Coarsely  breadcrusted  bombs  are 
relatively  dense,  with  well  defined,  dark-grey-to-black,  poorly-to-non- 
vesicular  rinds,  broad,  deep  surface  fractures  defining  large  polygons, 
and  grey-to-brown  vesicular  interiors  (Fig.  2e-f).  Finely  breadcrusted 
bombs  are  less  dense,  with  diffuse,  pale  vesicular  rinds,  finer  polygonal 
networks  of  narrower,  shallower  surface  fractures,  and  paler,  commonly 
flow-banded  interiors  (Fig.  2g-h).  Some  bombs  that  broke  during 
eruption  exhibit  two  generations  of  breadcrusting,  the  breakage  surface 
being  more  finely  breadcrusted  than  the  original,  outer  surface  of  the 
bomb.  Breakage  is  inferred  to  have  exposed  the  already  vesicular 
interior,  which  then  developed  a  second  generation  of  finer  breadcrust¬ 
ing.  Bombs  were  abraded  during  transport  in  the  pyroclastic  flows; 
most  lack  completely  preserved  breadcrust  surfaces  with  sharp  edges 
and  corners,  and  partial  rind  removal,  rounding  of  polygon  edges,  and 
abrasion  of  vesicular  interiors  are  common. 

Clasts  of  black,  essentially  nonvesicular  lava  resembling  the  glassy 
rinds  of  the  coarsely  breadcrust  bombs  are  interpreted  as  an  integral 
component  of  the  explosion-pyroclast  suite.  On  the  other  hand,  grey- 
to-brown  holocrystalline  lava  and  cinderblock  clasts  resembling 
typical  dome  rock  are  probably  derived  either  from  the  crater  walls 
or  from  earlier  block-and-ash  flow  deposits  traversed  by  the 
explosion  pyroclastic  flows. 


Table  1 

Results  of  textural  analysis  of  explosion  pyroclasts. 


BCP1  rind 

BCP43  rind 

BCP1  core 

BCP43  core 

AM029 

AM036 

PV3 

R2 

Sample  type 

Coarsely 
breadcrusted 
bomb  rind 

Finely 

breadcrusted 
bomb  rind 

Coarsely 

breadcrusted 

bomb  core 

Finely 

breadcrusted 

bomb  core 

Pyroclastic 
flow  pumice 
(dense) 

Pyroclastic 
flow  pumice 
(dense) 

Pyroclastic  flow 
pumice  (highly 
vesicular) 

Fallout  pumice 
(dense) 

Bulk  vesicularity  —  He-pycnometry  (vol.%) 

8±2 

12  ±  2 

33  ±2 

49  ±1 

58  ±1 

61  ±1 

76  ±1 

- 

Isolated  fraction  —  He-pycnometrya 

1.00 

- 

0.33 

0.13 

0.15 

0.09 

0.09 

- 

Bulk  vesicularity  —  image  analysis  (vol.%) 

9±2 

9  ±  2 

36±2 

48  ±2 

61  ±2 

57  ±2 

74  ±  2 

56  ±  2 

Isolated  fraction—  image  analysis3 

0.13 

0.37 

0.23 

0.46 

0.09 

0.07 

0.11 

0.17 

Crystallinity  vesicle-free  —  image  analysis  (vol.%) 

54 

59 

79 

54 

73 

66 

55 

69 

Glass  —  image  analysis  (vol.%) 

46 

41 

21 

46 

27 

34 

45 

31 

Modes  of  vesicle  populations.13  (pm) 

Small 

10 

14 

20 

20 

23 

16 

24 

13 

Intermediate 

- 

- 

- 

142 

82 

40 

135 

80 

Large 

1273 

455 

1221 

2299 

1879 

735 

2145 

1432 

Proportions  of  vesicle  populations.5 

Small 

15 

37 

62 

71 

40 

42 

38 

53 

(vol.%  of  bulk  vesicularity) 

Intermediate 

- 

- 

- 

9 

31 

44 

25 

34 

Large 

85 

63 

38 

20 

29 

12 

37 

13 

Nv  all  vesicles  (x  1014  m“ 3  of  glass)c 

Cheng  and  Lemlich  (1983) 

9.5  ±0.4 

26.2  ±1.1 

22.6  ±1.2 

48.0  ±2.8 

15.9  ±1.0 

31.1  ±3.7 

23.9  ±2.0 

66.8  ±  4.9 

Nv  connected  vesicles  (x  1014  m“ 3  of  glass)c 

0.0  ±0.0 

0.0  ±0.0 

17.9  ±1.0 

18.6  ±1.2 

7.5  ±0.6 

22.9  ±3.4 

17.8  ±1.4 

48.8  ±4.1 

Nv  isolated  vesicles  (x  1014  m“ 3  of  glass)0 

9.5  ±0.4 

26.2  ±1.1 

4.6  ±0.2 

29.4  ±1.6 

8.3  ±  0.4 

8.2  ±0.3 

6.1  ±0.7 

18.0  ±0.9 

Nv  crystals  (x  1014  m“ 3  of  glass)0 

139±4 

162  ±4 

357  ±9 

137  ±  4 

177  ±.5 

170  ±  4 

131  ±  3 

357  ±8 

Nv  oxides  (x  1014  m“ 3  of  glass)0 

83  ±2 

101  ±3 

174±  5 

89  ±3 

77  ±2 

82  ±2 

97  ±2 

163  ±4 

Nv  all  vesicles  (x  1014  m“ 3  of  glass)0 

Sahagian  and  Proussevitch  (1998) 

9.9  ±0.4 

26.8  ±1.1 

21.7  ±  1.3 

43.4  ±3.0 

15.9  ±  1.1 

27.8  ±4.3 

25.7  ±2.5 

57.5  ±5.1 

Nv  connected  vesicles  (x  1014  m“ 3  of  glass)0 

0.0  ±0.0 

0.0  ±0.0 

17.1  ±1.1 

18.0  ±1.3 

7.4  ±0.6 

20.2  ±4.0 

19.1  ±1.7 

41.2  ±4.2 

Nv  isolated  vesicles  (x  1014  m“ 3  of  glass)0 

9.9  ±0.4 

26.8  ±1.1 

4.6  ±0.2 

25.4  ±1.7 

8.5  ±0.5 

7.6  ±0.3 

6.6  ±0.9 

16.3  ±0.9 

Nv  crystals  (x  1014  m“ 3  of  glass)0 

147  ±3 

163  ±3 

413  ±  6 

142  ±3 

177  ±.4 

182  ±  3 

143  ±2 

407  ±6 

Nv  oxides  (x  1014  m“ 3  of  glass)0 

96  ±1 

113  ±2 

208  ±3 

101  ±  2 

78  ±1 

95  ±1 

112  ±  1 

198  ±  3 

Vesicles  number  of  each  population. 

Small 

9.5xl014 

26.2  xlO14 

22.6xl014 

48.0x  1014 

15.7  x  1014 

30.8  xlO14 

23.9xl014 

66.6xl014 

(m~3  of  melt) 

Intermediate 

- 

- 

- 

1.1  x  1011 

1.4xl013 

3.8  xlO13 

0.4xl013 

1.2xl013 

Large 

1.5xl010 

2.5  xlO10 

8.3  x  1010 

3.3  xlO8 

2.1  xlO9 

9.9  xlO9 

9.0  xlO9 

2.3  xlO9 

Fractal  dimension 

-1.9 

-2.3 

-2.7 

-2.8 

-3.1 

-3.4 

-2.9 

-3.3 

H20  in  residual  glass  (wt.%)d 

0.83  ±0.13 

- 

1.10  ±0.34 

- 

0.45  ±0.31 

0.74  ±0.33 

1.32  ±0.14 

0.75  ±  0.13e 

Quench  pressures  Psat  (MPa)  from  H20  content01 

6.1  ±1.8 

- 

10.3  ±5.9 

- 

1.9  ±2.5 

4.9  ±4.1 

14.5  ±2.8 

5.0±1.5e 

Pre-explosive  vesicularity  (vol.%)  estimated  by 
-  Vesicle  size  populations1 

<1 

<1 

<1 

<1 

28 

11 

47 

12 

-  Recorded  strain  by  amphiboles1 

- 

- 

- 

- 

35 

27 

57 

27 

3  Isolated  fraction  =  isolated  vesicularity /total  vesicularity. 
b  Determined  using  the  KWare  SFT  software. 

c  The  errors  assume  measurements  of  vesicle  numbers  to  within  ±  5%  on  individual  images. 
d  From  Burgisser  et  al.  (in  press). 

e  Unpublished,  from  S.  Poussineau  (personal  communication). 
f  See  text  for  explanation. 
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gradient  in  vesicularity,  line  marks  the  outer  surface  of  the  clast,  d)  Curviplanar  tears  and  cracks  subparallel  to  clast  margins  in  a  dense  pumice,  e-f)  Exterior  and  cross  section  of  a 
coarsely  breadcrust-bomb,  g-h)  Exterior  and  cross  section  of  a  finely  breadcrusted  bomb. 
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5.  Pyroclast  vesicularities 

Vesicularities  of  texturally  homogenous  pumice  lapilli  and  blocks 
range  from  24  to  79vol.%  (Fig.  3)  and  correlate  with  colour,  being 
lowest  in  darker  pumices  and  higher  in  paler  ones.  The  fraction  of 
isolated  vesicles  (isolated  divided  by  total  vesicularity)  is  universally 
low  (<0.25,  with  85%<0.1 ).  Flow  pumices  cover  the  entire  vesicularity 
range  and  have  isolated  fractions  of  0-0.13,  whereas  fallout  pumices 
have  vesicularities  of 43-72  vol.%  and  isolated  fractions  of  0.04-0.14,  a 
single  sample  having  0.25  (Fig.  3).  No  variation  of  either  vesicularity 
or  isolated  vesicle  fraction  with  clast  size  is  observed. 

Vesicularity  profiles  across  eight  >30  cm  pumice  blocks  are  shown 
in  Fig.  4.  Four  of  these  appeared  homogeneous  in  the  field,  and  four 
had  visually  obvious  radial  gradients  in  vesicle  size.  The  four 
homogeneous  blocks  (SHV4-1 2-1 3-22)  lack  significant  gradients  in 
vesicularity  from  core  to  rim,  as  anticipated  from  inspection.  The  four 
vesicle-size-graded  blocks  (SHV2-14-23-25),  on  the  other  hand, 
exhibit  vesicularity  gradients,  but  these  vary  from  sample  to  sample 
and  no  systematic  decrease  in  vesicularity  from  core  to  rim  is  evident. 
The  coarse  interiors  of  these  pumices  are  no  more  vesicular  than  the 
more  finely  vesicular  rims.  Textural  coarsening  in  the  interiors 
therefore  took  place  without  inflation,  as  consistent  with  the  absence 
of  surface  breadcrusting. 

Breadcrust-bombs  differ  from  pumices  in  that  (1)  their  vesicular¬ 
ity  range  (20-66  vol.%  most  lying  between  35  and  55%,  Fig.  3)  is 
smaller,  and  highly  vesicular  (>66%)  samples  are  not  observed;  and 
(2)  the  fraction  of  isolated  pores  (0.05-0.33,  >80%  being  0.1 -0.2) 
is  higher  than  in  pumices  of  similar  vesicularity.  Bomb  rinds  contain 
1-25  vol.%  vesicles,  most  of  which  are  isolated.  Rind  and  interior 
vesicularities  are  broadly  correlated  (Fig.  5).  Coarsely  breadcrusted 
bombs  have  the  lowest  vesicularities,  both  in  rinds  and  interiors,  and 
finely  breadcrusted  bombs  are  more  vesicular.  It  is  the  existence  of 
vesicular  rinds  on  finely  breadcrusted  bombs  that  gives  these  bombs 
their  pale  colours  and  make  distinction  between  rind  and  interior  less 
clear  than  in  the  coarsely  breadcrusted  bombs. 

Full  tables  of  vesicularity  data  are  provided  as  Supplementary 
electronic  material. 

6.  Microscopic  vesicle  textures 

The  pyroclasts  contain  vesicles  with  a  broad  range  of  sizes  set  in 
microlite-bearing  groundmass.  In  this  section  we  focus  on  vesicles  less 
than  a  few  mm  in  diameter  present  in  hand  specimens,  and  distinguish 
three  populations:  small  (less  than  a  few  tens  of  pm),  intermediate  (few 
tens  to  a  few  hundreds  of  pm)  and  large  (few  hundreds  of  pm  to  a 
few  mm).  It  is  shown  later  that  these  three  populations  also  have  genetic 
significance. 


Fig.  3.  Plot  of  connected  versus  total  vesicularity  for  all  the  samples  of  this  study.  Lava 
samples  from  dome  collapses  at  Soufriere  Hills  Volcano  are  also  shown  (Formenti  and 
Druitt,  2003).  As  connected  vesicularities  could  not  be  determined  for  breadcrust 
bombs  rinds,  we  just  show  the  range  of  bulk  vesicularities  obtained  (thick  black  line). 


Fig.  4.  Vesicularity  as  a  function  of  relative  position  inside  large  pumices  (>30  cm). 
Filled  diamonds  with  solid  lines  are  those  pumices  that  were  judged  in  the  field  to  be 
texturally  homogeneous:  squares  with  dashed  lines  are  those  that  had  larger  vesicles  in 
the  interior  than  in  the  rind. 


Vesicle  textures  in  fallout  and  flow  pumices  are  very  similar  and  are 
described  together.  The  large  vesicles  form  interconnected  networks 
with  curved,  scalloped  walls  indicative  of  coalescence.  Large  vesicles  in 
the  more  vesicular  pumices  are  quasi-spherical  to  elliptical  in  shape. 
Those  in  dense  pumices  commonly  have  more  ragged,  fissure-like 
shapes,  suggesting  that  perhaps  they  already  existed  prior  to  explosion. 
About  15%  of  the  large  vesicles  are  angular  voids  associated  with 
fractured  amphibole  phenocrysts  (Fig.  6a).  Intermediate-sized  vesicles 
in  all  pumices  have  variably  rounded  to  ragged  shapes  and,  like  the  large 
ones,  form  interconnected  networks  in  three  dimensions.  In  contrast, 
small  vesicles  are  commonly  spherical  and  many  are  isolated;  they 
either  form  a  ‘matrix’  in  which  the  intermediate  vesicles  are  dispersed 
(Fig.  6b),  or  are  situated  in  the  walls  separating  the  latter.  In  some 
samples  the  smallest  isolated  vesicles  form  sub-spherical  clusters 
several  tens  of  microns  in  diameter  that  protrude  with  bulbous, 
cauliform  shapes  into  larger  vesicles  (Fig.  6c;  Formenti  and  Druitt, 
2003).  There  is  textural  evidence  that  many  vesicles  of  intermediate  size 
formed  by  coalescence  of  the  small  vesicles  (rather  than  pre-existing 
them),  the  process  commonly  being  preserved  quenched  in  progress 
(Fig.  6d).  The  sizes  of  some  intermediate  vesicles  appear  to  be  inherited 
from  the  clusters  of  small  vesicles  when  the  latter  coalesced  while 
preserving  the  overall  sub-spherical  form  of  the  cluster. 

Vesicles  in  pumices  are  commonly  observed  in  spatial  association 
with  phenocrysts.  Large,  angular  voids  are  associated  with  fractured 
amphiboles,  and  have  two  endmember  types:  (1)  voids  in  amphiboles 
boudinaged  uniaxially  in  the  plane  of  flow  foliation,  with  well  defined 
length-perpendicular  fractures  (Fig.  7);  (2)  voids  in  amphiboles  that  are 
fractured  both  perpendicular  and  parallel  to  length,  and  the  fragments 
dispersed  around  the  vesicle  margins  in  a  manner  suggestive  of  more 
isotropic  expansion.  In  both  types,  crystal  fragments  are  commonly 


Fig.  5.  Relationships  between  rind  and  interior  vesicularities  of  breadcrust  bombs, 
including  both  coarsely  and  finely  breadcrusted  types. 
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Fig.  6.  SEM  images  of  broken  surfaces  (a-d)  and  thin  sections  (e-f)  of  pumices,  a)  Angular  void  in  a  fractured  amphibole  phenocryst,  the  fragments  being  connected  by  thin  glass 
fibres  (white  arrows),  b)  Visual  evidence  for  three  different  size  populations  (large,  intermediate  and  small)  of  vesicles  in  pumices,  c)  Cauliform-shaped  clusters  of  small  vesicles 
protruding  into  intermediate  ones,  d)  Evidence  for  coalescence  of  small  vesicles  to  form  intermediate-sized  ones,  e)  Microphenocryst  of  plagioclase  surrounded  by  radiating, 
elongated  vesicles,  f)  Wrinkling  of  vesicle  wall  indicative  of  the  onset  of  rupture  (white  arrow). 


connected  by  thin,  delicate  threads  of  glass  generated  either  by  the 
bursting  of  melt  inclusions,  or  by  the  pulling-out  of  thin,  pre-existing 
melt  films  in  incipient  cracks.  A  single  type  of  amphibole-associated  void 
is  commonly  dominant  within  a  given  pumice  block.  Type  1  is  observed 
in  -45%  of  pumices  and  type  2  in  -35%,  the  remaining  -  20%  of  pumices 
lacking  voids  associated  with  amphibole.  Another  common  texture 
involves  radial  arrangements  of  stretched  vesicles  around  phenocrysts 
of  plagioclase  or  amphibole  (Fig.  6e).  This  is  attributed  to  expansion  of  a 
magmatic  foam  around  a  rigid  crystal;  it  cannot  be  due  to  heterogeneous 
bubble  nucleation  because  in  each  case  the  vesicles  are  separated  from 
the  crystal  by  a  thin  glass  film,  showing  that  the  crystal  was  not  wetted 
by  gas.  Only  in  the  case  of  titanomagnetite  is  it  common  to  see  vesicles  in 
direct  contact  with  crystals  without  intervening  glass,  suggesting  that 
titanomagnetite  provided  nucleation  sites  for  bubbles  (Fig.  8). 

There  is  abundant  evidence  that  bubble  coalescence  was  ongoing  at 
all  scales  larger  than  a  few  pm  at  the  time  of  sample  quench:  ovoid,  neck¬ 
like  connections  with  partially  retracted  walls  between  neighbouring 
vesicles  (Fig.  6d),  wrinkling  of  thin  vesicle  walls  (Fig.  6f),  the  occurrence 
of  thin  glass  fibres,  and  the  interconnection  of  all  but  a  fraction  of  the 
smallest  vesicles.  Minimum  observed  vesicle  wall  thicknesses  are  <  1  pm. 

Breadcrust  bomb  rinds  contain  small,  mostly  isolated,  vesicles  that 
are  irregularly  distributed,  being  most  abundant  near  rind-penetrating 


surface  fractures  and  around  phenocrysts  (Fig.  9a,  c).  Areas  of  vesicle- 
free  groundmass  occur  in  the  rinds  of  coarsely  breadcrusted  bombs,  but 
not  in  those  of  the  finely  breadcrusted  bombs.  The  lower  limit  of  the  rind 
is  commonly  marked  by  string-like  networks  of  small  vesicles,  which 
then  merge  to  form  the  more  uniformly  distributed  vesicle  population  of 
the  interior.  The  interiors  of  all  bombs  contain  distinct  large  and  small 
vesicle  populations.  Large  vesicles  are  invariably  associated  with 
fractured  amphiboles,  like  those  in  the  pumices.  However,  well 
developed  uniaxial  boudinage  is  never  observed  in  breadcrust  bombs, 
and  the  voids  are  mostly  of  the  more  isotropic  type  2.  Small  vesicles  are 
uniformly  distributed  throughout  the  bomb  interiors  (Fig.  9b,  d);  they 
are  mostly  isolated,  with  quasi-spherical  forms,  and  commonly  occur  in 
strings  and  clusters  around  crystals  and  large  vesicles.  Evidence  for 
vesicle  coalescence  is  abundant  in  bomb  interiors,  although  less  so  than 
in  pumices. 

7.  Size  distributions  of  vesicles  and  crystals 

The  six  samples  chosen  for  analysis  of  vesicle  and  crystal  size 
distributions  were  a  coarsely  breadcrusted  bomb  (BCP1),  a  finely 
breadcrusted  bomb  (BCP43),  three  pyroclastic-flow  pumices  (AM029, 
AM036  and  PV3),  and  a  fallout  pumice  (R2).  Separate  measurements 
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Fig.  7.  Uniaxially  boudinaged  amphibole  phenocryst  in  pumice  indicative  of  significant 
syn-explosive  vesiculation  prior  to  magma  fragmentation. 


were  made  on  the  rinds  and  interiors  of  the  two  breadcrust  bombs. 
Groundmass  water  contents  of  five  of  the  samples  were  analysed  by 
Burgisser  et  al.  (in  press).  The  data  are  listed  in  Table  1. 

Glass-referenced  vesicle-volume  distributions  are  shown  in  Fig.  10, 
with  connected  and  isolated  vesicles  distinguished.  The  limits  be¬ 
tween  different  image  magnifications  are  also  shown.  No  major  artefacts 
related  to  changes  in  magnification  are  apparent;  some  small  dis¬ 
continuities  in  the  distributions  of  AM029  and  PV3  may  be  artefacts,  but 
no  conclusions  are  based  on  these.  Total  vesicularities  calculated  by 
image  analysis  are  within  4  vol.%  of  pycnometry  measurements  on  the 
same  clasts  (Table  1). 

Vesicle-volume  distributions  of  the  breadcrust  bombs  are  mark¬ 
edly  bimodal.  The  large  population  is  accounted  for  entirely  by 
amphibole-associated  voids.  The  small  vesicles  form  a  narrow,  well 
defined  peak;  those  in  the  rinds  are  entirely  isolated  and  form  a  mode 
at  <10  pm,  whereas  in  the  bomb  interiors  a  significant  fraction  are 
interconnected  and  the  size  mode  is  larger  (-20  pm).  The  interior  of 
the  finely  breadcrusted  bomb  contains  a  minor  population  of 
intermediate-sized  vesicles.  The  distributions  of  all  four  pumices  are 
also  roughly  bimodal,  with  a  size  gap  at  -  200-500  pm.  The  finer  mode 
is,  however,  much  broader  than  in  the  breadcrust  bombs,  consisting  of 
comparable  proportions  of  small  and  intermediate-sized  vesicles.  The 
small  population  resembles  that  in  the  breadcrust  bombs  and  contains 
all  the  isolated  vesicles  present.  The  intermediate  population  forms  a 
coarse-tail  shoulder  on  the  small  population.  In  pumices  PV3  and  R2, 


Fig.  8.  Textural  evidence  for  heterogeneous  nucleation  on  titanomagnetite.  The  crystal 
appears  to  lack  a  glass  rind,  indicating  that  the  gas  wetted  the  solid.  This  is  not  the  case  for 
amphibole  (Fig.  6a)  or  plagioclase  (Fig.  6e),  which  always  have  thin  glass  films  separating 
crystal  from  vesicle.  The  wetting  angle  of  this  magnetite  crystal  is  -130°,  consistent  with 
the  high  wetting  angle  (158.5°)  proposed  by  Hurwitz  and  Navon  (1994)  for  magnetite. 


the  small  and  intermediate  populations  form  distinct  sub-peaks  at 
10-20  and  65-80  pm,  respectively. 

The  proportions  of  the  three  populations  were  estimated  from 
Gaussian  best  fits  using  the  KWare  SFT  software.  There  is  no  physical 
basis  for  this;  it  is  simply  better  than  eyeballing.  Grouping  all  samples, 
the  modes  of  the  three  populations  are  10-24  pm  (small),  40-142  pm 
(intermediate)  and  455-2299  pm  (large)  (Table  1).  The  small 
population  constitutes  15-37%  of  the  vesicularity  of  the  breadcrust 
bomb  rinds,  62-71%  of  that  of  the  breadcrust  bomb  interiors  and  38- 
53%  of  that  of  the  pumices.  The  intermediate  and  large  populations 
account  for  25-44%  and  12-37%  of  pumice  vesicularity,  respectively. 

Vesicle  size  distributions  can  also  be  plotted  using  vesicle  population 
number  density,  which  is  equal  to  the  slope  of  the  cumulative  number 
distribution  (Fig.  11).  Number  distributions  are  dominated  by  the  smallest 
vesicles;  in  our  samples  more  than  99%  of  the  vesicularity  in  number  terms 
consists  of  vesicles  of  the  small  population.  The  distributions  are  curved  for 
all  samples  and  at  all  size  scales,  emphasizing  the  importance  of 
coalescence  of  vesicles  larger  than  a  few  pm.  If  number  distributions  are 
redrawn  cumulatively  (Fig.  12),  then  all  samples  define  approximately 
linear  trends  at  >30  pm.  Linear  trends  on  such  plots  depict  power-law 
distributions,  the  slope  being  a  fractal  dimension  (Gaonac’h  et  al.,  1996a,b; 
Blower  et  al.,  2001).  The  data  give  slopes  ranging  from  —1.9  to  —3.4, 
being  lowest  in  pumices  (  —  2.9  to  —3.4),  intermediate  in  breadcrust- 
bomb  interiors  (  —  2.7  to  —2.8)  and  highest  in  breadcrust-bomb  rinds 
( —  1 .9  to  —  2.3;  Table  1 ).  Vesicles  smaller  than  -30  pm  fall  off  the  power- 
law  trends,  defining  a  different  tendency.  Total  vesicle  number  densities 
range  from  10  to  67  x  1014  m~3  (method  of  Cheng  and  Lemlich  1993). 
Values  for  isolated  vesicles  alone  range  from  10  to  29xl014m~3.  No 
systematic  correlation  with  pyroclast  vesicularity  is  observed. 

A  vesicle  circularity  index  C=4nA/P2  (A  is  surface  area  and  P  is 
perimeter)  is  plotted  on  Fig.  10  as  a  function  of  vesicle  size,  with  separate 
curves  for  isolated  and  connected  vesicles.  A  spherical  vesicle  would 
have  a  C^l,  whereas  a  chain  of  newly  coalesced  vesicles  would 
approach  C^0.  Isolated  vesicles  are  systematically  more  circular  than 
connected  ones  of  similar  size.  Despite  some  artefacts  related  to  image 
magnification  limits  (e.g.,  -40  and  -250  pm  in  AM029),  C  values  of 
connected  vesicles  in  general  decrease  with  increasing  vesicle  size 
throughout  the  intermediate  and  large  populations,  consistent  with 
their  more  ragged  shapes  inherited  from  coalescence,  the  number  of 
coalescence  events  increasing  with  vesicle  size. 

Crystal  volume  distributions  for  the  six  pyroclasts  (Fig.  13)  are 
presented  for  comparison  with  the  vesicle  distributions.  Two  main 
populations  can  be  recognised;  phenocrysts  larger  than  a  few 
hundred  pm,  and  microphenocrysts  and  microlites  of  -300  to  1  pm, 
a  bimodality  being  present  in  three  samples  (AM029,  AM036  and  R2). 
Total  crystal  contents,  calculated  vesicle-free,  are  lowest  (55  vol.%)  in 
the  most  vesicular  pumice  (PV3),  intermediate  in  denser  pumices 
(66-73  vol.%;  AM029,  AM036  and  R2),  and  highest  in  the  breadcrust 
bombs  (54-79  vol.%)  (Table  1).  This  is  due  primarily  to  a  greater 
abundance  of  microphenocrysts  and  microlites  in  the  denser 
pyroclasts  (Fig.  13),  as  observed  by  other  authors  (Murphy  et  al., 
1998;  Clarke  et  al.,  2007).  Total  glass-referenced  crystal  number 
densities  range  from  131  to  357xl014  m~3,  values  for  Fe-Ti  oxides 
(mainly  titanomagnetite)  alone  being  77  to  174xl014  m-3.  No  sys¬ 
tematic  variations  of  crystal  number  density  with  pyroclast  type,  or 
between  rinds  and  cores  of  breadcrust  bombs,  are  observed.  There  is 
no  correlation  between  crystal  and  vesicle  number  densities. 

8.  Discussion 

8.1.  Pre-explosion  conduit  stratigraphy 

Previous  studies  enable  us  to  reconstruct  conduit  stratigraphy 
immediately  prior  to  an  average  1997  explosion.  Approximate  conduit 
pressures  for  the  explosion  pyroclasts  have  been  estimated  using  two 
independent  methods,  microlite  abundance  and  groundmass-glass 
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Fig.  9.  Backscattered  SEM  images  of  breadcrust  bombs.  All  the  images  are  at  the  same  scale.  Labels  are:  (A)  amphibole,  (PI)  plagioclase,  (0)  oxide,  ( Vc)  connected  vesicle,  (Vi)  isolated 
vesicle,  a)  Coarsely  breadcrusted  bomb  rind  containing  rare  small  isolated  vesicles.  Note  that  the  bomb  rinds  are  typically  very  heterogeneous,  and  vesicles  are  not  present 
everywhere.  The  occurrence  of  areas  of  rind  free  of  vesicles  shows  that  the  coarsely  breadcrusted  bombs  inflated  from  essentially  non-vesicular  magma,  b)  Vesicular  interior  of 
coarsely  breadcrusted  bomb  core  containing  both  isolated  (non-impregnated:  grey)  and  connected  (impregnated:  black)  vesicles,  c)  Finely  breadcrusted  bomb  rind  with  patchily 
distributed  small  vesicles;  note  the  higher  vesicle  numbers  than  in  the  rind  of  the  coarsely  breadcrusted  bomb  (a),  d)  Finely  breadcrusted  bomb  core  containing  the  same  two  types 
of  vesicles  as  the  coarsely  breadcrusted  bomb  (b),  but  in  greater  abundance. 


water  content,  both  of  which  yield  similar  results  (Clarke  et  al.,  2007; 
Burgisser  et  al,  in  press).  Glassy  lava,  breadcrust  bombs  and  poorly 
vesicular  (<40  vol.%)  pumices  yield  PHi 0  values  of  0-10  MPa,  whereas 
more  vesicular  pumices  (>40  vol.%)  yield  up  to  55  MPa.  The  data  are 
consistent  with  the  existence  of  a  dominant,  gas-rich  zone  (highly 
vesicular  pumices)  that  passed  upwards  into  a  thin,  largely  degassed 
plug  (lava  fragments  and  breadcrust  bombs)  through  a  partly  degassed 
region  (dense  pumices)  (Clarke  et  al.,  2007;  Burgisser  et  al.  in  press). 
This  interpretation  ignores  any  possible  horizontal  variations  in  the 
conduit  (Kennedy  et  al.,  2005).  Variations  in  pyroclast  proportions 
between  different  pyroclastic  flow  lobes  are  attributed  to  withdrawal 
from  different  levels  in  the  conduit,  lobes  richer  in  dense  pumice  being 
derived  from  higher  levels  of  the  conduit  and  those  richer  in  more 
vesicular  pumice  coming  from  deeper. 

Crystals  (phenocrysts,  micro-phenocrysts  and  microlites)  in  the 
pyroclasts  are  interpreted  as  having  grown  prior  to  explosion,  either 
at  depth  or  during  slow  viscous  ascent  between  explosions  (Clarke 
et  al.,  2007).  This  is  supported  by  the  lack  of  correlation  between 
the  number  densities  of  crystals  and  vesicles  in  our  samples,  as  well  as 
by  the  presence  of  high  crystal  number  densities  in  the  (rapidly 
quenched)  rinds  of  coarsely  breadcrusted  bombs  (Table  1). 

8.2.  Origins  of  the  vesicle  populations 

A  vesicle  is  here  termed  ‘syn-explosive’  if  it  nucleated  during  and/or 
immediately  following  the  explosion  that  ejected  its  host  pyroclast.  It  is 
‘pre-explosive’  if  it  already  existed  prior  to  the  explosion  in  question, 
having  nucleated  either  during  previous  explosions  and/or  during  the 
3-63  hours  of  slow  viscous  ascent  between  explosions. 

Vesicles  of  the  small  population  (less  than  a  few  tens  of  pm)  are 
present  in  all  pyroclasts,  and  are  interpreted  as  having  nucleated  and 
grown  syn-explosively.  The  principal  argument  is  that  they  dom¬ 


inate  the  interiors  of  coarsely  breadcrusted  bombs,  in  which  a  syn- 
explosive  origin  is  unambiguous.  Coarsely  breadcrusted  bombs 
inflated  syn-explosively  from  essentially  non-vesicular  initial  states, 
as  shown  by  the  common  presence  in  the  rinds  of  large  areas  of 
vesicle-free  groundmass  (e.g.,  Walker,  1982;  Wright  et  al.,  2007).  All 
vesicles  in  the  bomb  interiors  must  therefore  be  syn-explosive  and 
post-date  fragmentation.  The  small  size,  high  sphericity  and  high 
isolated  fraction  of  these  vesicles  is  consistent  with  nucleation  and 
rapid  growth  shortly  before  quenching  (e.g.  Mourtada-Bonnefoi  and 
Laporte,  2002;  Mangan  et  al.,  2004;  Gardner  and  Denis,  2004).  The 
similarity  between  small  vesicles  in  the  pumices  and  finely  bread¬ 
crusted  bombs  with  those  in  the  coarsely  breadcrusted  bombs  leads 
us  to  interpret  all  vesicles  of  this  population  as  syn-explosive. 

Vesicles  of  the  intermediate  population  occur  in  all  pumices  as  a 
coarse  tail  to  the  small  population;  they  are  also  present  in  minor 
quantities  in  the  interiors  of  finely  breadcrusted  bombs,  but  are  absent 
from  coarsely  breadcrusted  bombs.  Since  many  are  observed  to  have 
formed  by  coalescence  of  the  small-population  vesicles,  they  too  must 
be  largely  syn-explosive  in  origin.  Some  appear  to  have  inherited  their 
sizes  from  internally  coalesced  vesicle  clusters,  perhaps  explaining  the 
distinct  small  and  intermediate  size  peaks  in  some  pumices  (Fig.  10). 

Large-population  vesicles  are  interpreted  as  having  a  mixed  origin. 
Those  in  the  breadcrust  bombs  are  invariably  amphibole-associated 
voids  that  must  be  syn-explosive  because  (1)  they  are  present  in  the 
interiors  of  coarsely  breadcrusted  bombs,  but  not  in  the  rinds,  (2)  they 
are  commonly  traversed  by  delicate  glass  filaments  indicative  of  a 
very  late-stage  origin.  The  -15%  of  large  vesicles  associated  with 
boudinaged  amphiboles  in  pumices  are  also  inferred  to  be  syn- 
explosive,  for  the  same  reasons.  The  remaining  -85%  of  large  vesicles 
in  pumices  are  not  associated  with  phenocrysts,  and  may  be  pre¬ 
explosive  in  origin,  although  possibly  enlarged  by  syn-explosive 
decompression.  The  pre-explosive  origin  of  these  large  vesicles  is 
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Fig.  11.  Vesicle  size  distributions  for  the  eight  pyroclasts  analysed  in  detail. 


Vesicle  Equivalent  Diameter  D  Qim) 

Fig.  12.  Cumulative  number  densities  ( NV>D )  for  the  eight  pyroclasts  analysed  in  detail. 

supported  by  the  abundance  in  some  dense  pumices  of  fissure-like 
networks  of  large  vesicles  that  appear  to  be  too  well  developed  to 
have  formed  during  explosion,  and  which  are  tentatively  interpreted 
as  having  acted  as  outgassing  pathways  prior  to  explosion. 

Some  cm-sized  ductile  tears  and  cracks  parallel  to  flow  banding  or 
clast  margins  commonly  observed  in  large  pumice  blocks  are 
envisaged  to  have  formed  (syn-explosively)  upon  ballistic  impact  or 
during  transport  in  pyroclastic  flows. 

Ignoring  second-order  modes,  the  vesicle-size  distributions  of  all 
pumices  and  breadcrust-bomb  interiors  are  essentially  bimodal 
(Fig.  10),  with  a  gap  from  a  hundred  to  a  few  hundred  pm.  Vesicle 
bimodality  is  observed  commonly  in  pyroclasts  from  explosive 
eruptions  (Witham  and  Sparks,  1986;  Klug  et  al.,  2002).  In  the 
breadcrust  bombs  of  the  present  study,  bimodality  is  explained  by 
syn-explosive  vesiculation  at  two  scales:  (1 )  nucleation  and  growth  of 


small  bubbles  and  (2)  formation  of  large  voids,  the  sizes  of  which  are 
imposed  by  those  of  amphibole  phenocrysts.  Bimodality  in  the 
pumices  is  tentatively  explained  by  combining  (1)  vesicles  formed 
by  syn-explosive  nucleation,  growth  and  coalescence  (small  and 
intermediate  populations)  with  (2)  a  large  population  composed  of 
syn-explosive,  phenocryst-related  voids  and  pre-explosive  bubbles. 

Interestingly,  the  size  gaps  in  the  pumice  vesicle-size  distributions 
correspond  approximately  with  gaps  in  corresponding  crystal  size 
distributions.  Crystals  in  all  four  pumices  have  size  gaps  between 
microlite/microphenocryst  and  phenocryst  populations  that  coincide 
approximately  with  those  of  vesicles.  This  is  a  real  feature,  not  a 
measurement  artefact.  Since  the  short  explosion  timescale  does  not 
permit  significant  syn-explosive  crystal  growth,  any  interrelationship 
must  involve  length-scale  imposition  on  growing  bubbles  by  crystals 
that  were  already  present.  Crystals,  being  rigid  objects,  will  influence 
the  sizes,  shapes  and  distributions  of  nearby  bubbles  (e.g.,  Fig.  6e). 
Voids  created  by  crystal  fracture  will  have  length  scales  imposed  on 
them  by  the  crystal  (e.g.,  Fig.  6a). 

8.3.  Amounts  of  pre- explosive  and  syn-explosive  vesiculation 

Syn-explosive  vesicularity  includes  the  small  and  intermediate 
populations,  all  the  large  population  in  breadcrust  bombs,  and  at  least 
~  1 5%  of  the  large  population  in  pumices.  The  breadcrust  bombs  inflated 
from  essentially  non-vesicular  magma,  whereas  the  pumices  had  pre¬ 
explosive  porosities  corresponding  to  no  more  than  -85%  of  the  large 
population,  quantified  as  follows.  Define  the  following  component 
fractions:  Xlv  (large  vesicles),  Xsiv  (small  +  intermediate  vesicles),  Xxt 
(crystals)  andXgl  (glass).  Pre-explosive  pumice  vesicularity  is  given  by 
1 00  x  0.85Xlv /  ( 0.85Xlv  +  Xxt  +  Xgl ) .  This  yields  11-47  vol.%  for  the 
pumices  in  Table  1,  and  is  greatest  for  the  most  vesicular  one  (47%; 
PV3).  These  are  maximum  estimates,  because  the  pre-explosive  bubbles 
could  have  expanded  during  explosion  decompression. 

An  independent  estimate  is  provided  by  uniaxially  boudinaged 
(type  1)  amphibole  phenocrysts  in  pumice,  which  are  interpreted  as 
providing  a  measure  of  strain  during  syn-explosive  expansion  (Fig.  7). 
Although  the  strain  recorded  by  some  isotropically  fractured  (type  2) 
amphiboles  could  be  partly  due  to  gas  release  from  vesiculated  melt 
inclusions,  we  consider  it  unlikely  that  this  is  the  case  for  the  uniaxial 


Fig.  10.  Glass-referenced  (i.e.,  crystal-free)  vesicle-volume  distributions  (lower  part  of  each  graph)  and  circularities  (upper  part  of  each  graph)  as  functions  of  size.  Connected 
vesicles  are  shown  in  grey  and  isolated  ones  in  black.  The  vertical  dashed  lines  are  the  magnification  limits  of  the  analysed  images.  Circularities  of  isolated  vesicles  marked  as  dashed 
lines  and  connected  vesicles  as  solid  lines.  The  black  arrows  on  each  graph  show  the  mode  of  each  of  the  small  (S),  intermediate  (I)  and  large  (L)  vesicle  populations  (Table  1). 


226 


T.  Giachetti  et  al.  /  Journal  of  Volcanology  and  Geothermal  Research  193  (2010)  215-231 


8 


6- 


2- 


0 


Coarsely  breadcrusted  bomb 
rind  (BCPlr) 


.pd 


DJ 


i 


i 


rt 


100 


u 


OpD 


1000 


10000 


8i 

6 

4- 

?■ 


iii  i 

Finely  breadcrusted  bomb 
rind  (BCP43r) 


g 

o 

>  0 

c3 

£81 


-  I 


■pD 


ood 


10 


a 


100 


1000 


10000 


u 


6- 


2^ 

0 

8 

6-1 


pnrffl 


aU 


10 


00 


2- 


0 


-pa 


•  innnflh 


n 


m 


Dense  pumice 
(AM029) 


a 


n 


1000 


10000 


Light  pumice 
(PV3) 


100 


D 


8 


oy 1 

8 

6 


oDDOIM 


Coarsely  breadcrusted  bomb 
core  (BCPlc) 


10 


100 


0 


1000 


10000 


0J_pd 
8 

6 


nil 


i 


Finely  breadcrusted  bomb 
core  (BCP43c) 


10 


oDDo 


i 


100 


a 


1000 


10000 


2- 


oLpd 


oDQdK 


Dense  pumice 
(AM036) 


00 


lop 


onO 


L 


0 


1000 


10000 


1 000  1 0000 

Crystal  Equivalent  Diameter  (pm) 


Fig.  13.  Glass-referenced  crystal  volume  distributions. 


case.  The  apparent  strain  varies  widely  between  pyroclasts,  probably 
due  to  different  amounts  of  syn-explosive  vesiculation  as  well  as 
variable  degrees  of  slip  between  melt  and  solid  at  the  crystal  surface, 
but  the  largest  values  (void  length  divided  by  crystal  length)  are  about 
2/3.  Now  defining  Xpre  as  the  fraction  of  pre-explosion  vesicles,  then 
Xjv+Xsiv— Xpre  is  the  fraction  of  syn-explosive  vesicles  and  we  have 
(Xiv + Xsiv  - Xpre)/ (Xpre +Xxt+Xgi)^  2/3.  Solving  this  for  Xpre  and 
converting  to  vesicularity  yields  maximum  pre-explosive  porosities  of 
27-57%,  which  are  similar  to  the  previous  estimates  and  again 
greatest  for  the  most  vesicular  pumice  (57%).  Amphibole  boudinage  is 
therefore  consistent  with  the  small  and  intermediate  vesicle  popula¬ 
tions  in  pumices  being  syn-explosive  in  origin.  We  conclude  that  pre- 
exposive  porosity  in  the  deep,  gas-rich  zone  of  the  conduit  was  <55%. 

8.4.  Timing  of  syn-explosive  vesiculation 

Breadcrust  bombs  and  pumices  preserve  evidence  for  different 
vesiculation  histories.  Bubble  nucleation  and  growth  in  the  pumices 
was  essentially  confined  to  the  conduit,  based  on  the  following 


evidence.  First,  the  absence  of  vesicularity  gradients  and  breadcrust- 
ing  in  pumices  shows  that  inflation  was  complete  before  surface 
quench  in  contact  with  the  air.  Second,  flow  and  fall  pumices  have 
similar  vesicle  textures,  contents  and  size  distributions,  despite 
different  thermal  histories  once  exited  the  vent.  Third,  the  existence, 
in  -45%  of  the  pumices,  of  amphiboles  boudinaged  uniaxially  in 
the  plane  of  flow  foliation  shows  that  vesiculation  of  at  least  this 
percentage  of  pumices  was  well  advanced  prior  to  fragmentation 
in  order  for  the  boudinage  to  be  uniaxial.  More  complex  modes  of 
amphibole  breakage,  both  perpendicular  and  parallel  to  crystal  length, 
may  be  due  to  increasingly  isotropic  decompression  and  inflation  as 
fragmentation  generated  progressively  smaller,  more  equant  magma 
parcels.  Syn-explosive  vesiculation  in  pumices  commenced  prior  to 
fragmentation,  and  in  many  cases  continued  after  fragmentation,  but 
was  complete  prior  to  the  onset  of  clast  surface  cooling.  The  peak  phase 
of  each  explosion  typically  lasted  a  few  tens  of  seconds  (Druitt  et  al. 
2002).  Once  in  the  air,  clast  surfaces  cooled  quickly  (another  10-20  s), 
causing  rapid  disappearance  of  incandescence  on  night  video  footage. 
Vesiculation  was  therefore  complete  in  less  than  a  minute  after 
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explosion  onset,  even  for  the  deepest-derived  magma.  Only  late-stage 
impact-generated  tearing  and  fracturing  of  pumices  generated  further 
voidage. 

Being  derived  from  the  top  of  the  conduit,  initially  non-vesicular 
magma  of  the  coarsely  breadcrusted  bombs  was  fragmented  and 
entered  the  air  almost  immediately  after  explosion  onset.  Immediate 
cooling  quenched  the  bomb  surface,  while  the  hot  core  inflated.  Near- 
universal  abrasion  of  rinds  and  interiors  shows  that  bomb  inflation 
was  complete  well  before  the  pyroclastic  flows  ceased  movement 
-300  s  later.  Detailed  work  on  breadcrust  bombs  at  Pinchincha 
Volcano  also  documented  rapid  bomb  inflation  (Wright  et  al.  2007). 
The  histories  of  finely  breadcrusted  bombs  are  inferred  to  have  been 
intermediate  between  those  of  pumices  and  coarsely  breadcrusted 
bombs.  Being  derived  from  slightly  deeper  in  the  conduit,  there  was 
some  time  for  initial  bubble  nucleation  and  growth  throughout  the 
clast  before  the  surface  was  quenched  and  further  inflation  was 
limited  to  the  interior. 

Syn-explosive  vesiculation  was  thus  controlled  by  the  time  available 
between  decompression  onset  and  quench  (Fig.  14).  This  was  shortest 
for  the  breadcrust  bomb  rinds,  intermediate  for  the  bomb  interiors,  and 
longest  for  the  pumices.  Neo-formed  vesicles  preserved  in  the  rapidly 
quenched  rinds  of  breadcrust  bombs  are  small  and  isolated.  During 
inflation  of  the  more  slowly  cooled  interiors,  these  vesicles  grew 
(shifting  the  mode  of  the  isolated  ones)  and  began  to  coalesce.  Syn- 
explosive  coalescence  is  most  advanced  in  the  pumices  because  the  time 
available  between  decompression  onset  and  quench  was  longest. 

8.5.  Mechanisms  of  syn-explosive  vesiculation 

8.5.1.  Bubble  nucleation 
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Fig.  14.  Schematic  illustration  showing  the  development  of  vesicle  volume  distributions 
of  breadcrust-bomb  rinds,  breadcrust-bomb  cores  and  pumices.  Isolated  syn-explosive 
vesicles  in  black  and  connected  syn-explosive  vesicles  in  dark  grey.  Those  in  light  grey 
are  angular  voids  formed  by  syn-explosive  fracturing  of  amphibole  phenocrysts.  Large, 
pre-explosive  vesicles  in  pumices  are  shown  in  white.  Coalescence  to  form  the 
intermediate-sized  vesicles  (dark  grey)  is  most  advanced  in  the  pumices. 


Homogeneous  nucleation  in  water-rich  rhyolitic  melts  requires 
supersaturation  pressures  of  60-310  MPa  (Mourtada-Bonnefoi  and 
Laporte,1999,  2002;  Mangan  and  Sisson,  2000),  which  exceed  the 
decompressions  of  a  few  to  -55  MPa  experienced  by  the  1997 
pyroclasts  (Burgisser  et  al.,  in  press).  This  implies  that  syn-explosive 
nucleation  took  place  heterogeneously.  Of  the  crystalline  phases 
present  in  Soufriere  Hills  magma,  only  titanomagnetite  has  been 
shown  to  provide  effective  sites  for  bubble  nucleation  at  super¬ 
saturations  as  low  as  a  few  MPa  (Hurwitz  and  Navon,  1994;  Gardner 
et  al.,  1999,  2000;  Gardner  and  Denis,  2004;  Gardner,  2007;  Cluzel 
et  al.,  2008).  Moreover,  Fe-Ti  oxides  (mostly  titanomagnetite)  are 
present  in  the  explosion  pyroclasts  with  number  densities  1.4  to  9 
times  higher  than  those  of  vesicles,  showing  that  sufficient  nucleation 
sites  were  available  (Table  1).  Textures  observed  on  SEM  images  are 
consistent  with  bubble  nucleation  on  titanomagnetite  (Fig.  8). 

The  number  densities  of  small  and  intermediate  vesicles  (Table  1 ) 
permit  estimation  of  syn-explosive  decompression  rates  (Toramaru, 
2006;  Cluzel  et  al.,  2008).  Nucleation  generates  a  stationary  bubble 
number  density  Nstat  that  is  independent  of  saturation  pressure,  but 
dependent  on  the  mean  decompression  rate  dP/dt: 


dP 

dt 


V/s 

Plj3) 


where  D  is  water  diffusivity,  a  is  interfacial  surface  tension  between 
melt  and  vapour,  Psat  is  initial  saturation  pressure,  T  is  absolute 
temperature,  and  a  is  a  constant  (3.5  xlO14  in  SI  units).  For  heteroge¬ 
neous  nucleation,  a  is  replaced  by  an  effective  value  in  the  presence  of  the 
crystalline  phase.  Heterogeneous  bubble  number  densities  Nv  theoret¬ 
ically  vary  both  with  the  number  of  available  nucleation  sites  and 
decompression  rate,  with  two  cases  (Cluzel  et  al.,  2008):  (1)  the  density 
of  available  nucleation  sites  Nn  is  greater  than  Nstat ,  in  which  case 
Nv~Nstat  and  Nv  depends  only  on  dP/dt  through  Eq.  (1 );  (2)  Nn<Nstaf,  in 
which  case  Nv  «  Nn  (<Nstat),  and  the  equation  cannot  be  used.  Since  oxide 
number  densities  in  our  samples  are  much  higher  than  vesicle  number 
densities  (case  1),  we  can  use  Nv  to  estimate  dP/dt.  This  assumes  that 
sufficient  time  was  available  for  bubble  number  densities  to  reach 
stationary  values  prior  to  quench;  if  not,  dP/dt  is  under-estimated. 

We  calculate  syn-explosive  dP/dt  values  for  a  typical  breadcrust 
bomb  interior,  a  dense  pumice  and  a  highly  vesicular  pumice  using 
our  data  for  Nv  (Table  1),  the  interstitial  glass  water  contents  and  Psat 
estimations  of  Burgisser  et  al.  (in  press),  and  cr~ 0.025  Nm”1 
appropriate  for  nucleation  on  titanomagnetite  (Cluzel  et  al.,  2008). 
The  values  range  from  0.3  to  6.5  MPa  s~  \  No  correlation  with  pyroclast 
type  is  observed,  perhaps  due  to  errors  associated  with  estimates  of 
water  content  (through  D),  Nv  and  Psat  for  the  different  samples. 

8.5.2.  Bubble  growth 

Bubble  growth  is  driven  by  decompression  and  volatile  diffusion, 
and  resisted  by  melt  viscosity  (Sparks,  1978;  Thomas  et  al.,  1994; 
Gardner  et  al.,  1996;  Proussevitch  and  Sahagian,  1996;  Gardner  et  al., 
2000).  Growth  during  decompression  from  an  initial  pressure  P,  is 
governed  by  three  characteristic  timescales  (Lensky  et  al.,  2004):  a 
decompression  time  Tdec  —  Psat/{ dP/dt),  a  diffusive  time  rdl^=Pf/D/ 
and  a  viscous  relaxation  time  rvisc  =  4rii/Psat ,  where  P„  D,-  and  q,  are 
bubble  radius,  diffusion  coefficient,  and  melt  viscosity  at  the  initial 
stages  of  growth.  If  Tvisc>Tdec,  then  viscosity  is  high  enough  to 
significantly  inhibit  bubble  growth,  whereas  if  Tvisc<^TdeCl  then  bubbles 
can  expand  freely  to  ambient  pressure.  The  latter  regime  can  in  turn  be 
subdivided.  If  diffusion  is  fast  relative  to  decompression,  rdlj<<c  Tdec,  then 
degassing  takes  place  at  equilibrium  and  bubbles  grow  by  both  diffusion 
and  decompression-driven  expansion,  whereas  if  the  converse  is  true, 
diffusive  transfer  is  negligible  and  bubbles  grow  by  decompression- 
driven  expansion  alone. 
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The  three  timescales  are  estimated  for  the  explosion  pyroclasts  in 
Table  2.  We  also  included  a  typical  breadcrust  bomb  rind  at  a 
temperature  close  to  that  of  the  glass  transition  (-700  °C  in  silicic 
melts  with  l-3wt.%;  Dingwell,  1998).  As  discussed  above,  decom¬ 
pression  of  a  breadcrust  bomb  from  the  top  of  the  conduit  took  place 
over  -10  s,  whereas  pumices  from  deeper  in  the  conduit  took  a  few 
tens  of  s.  Consider  the  case  of  a  10-pm-diameter  syn-explosive  bubble. 
For  the  pumices  we  find  that  rdec>Tdiff^>Tviso  implying  that  melt 
viscosity  was  sufficiently  low  for  bubbles  to  expand  freely  to  near¬ 
ambient  pressure.  Diffusion  was  able  to  occur  on  the  timescale  of 
decompression,  but  was  insufficiently  fast  to  attain  equilibrium,  in 
agreement  with  high  remnant  water  contents  in  interstitial  glasses. 
For  the  breadcrust  bomb  interior  Tdec>Tdijj>Tvisc,  implying  again 
disequilibrium  outgassing,  but  also  that  viscosity  was  high  enough  at 
these  lower  water  contents  to  now  partially  limit  bubble  growth. 
Finally,  cooling  of  the  breadcrust  bomb  rind  has  the  combined  effect  of 
significantly  raising  melt  viscosity  and  slowing  diffusion,  such  that 
TVisc^Tdiff>Tdec  and  any  bubbles  that  nucleated  were  unable  to  grow 
on  the  timescale  available  (i.e.,  they  were  viscosity-quenched; 
Thomas  et  al.,  1994),  again  in  agreement  with  observations. 

8.5.3.  Bubble  coalescence 

Textural  evidence  for  coalescence  is  found  in  all  vesicle  size  classes; 
only  the  small  population  contains  isolated  vesicles.  Bubble-wall  thinning 
and  rupture  in  an  expanding  foam  is  driven  by  bubble  growth  and 
pressure  gradients  (Navon  and  Lyakhovsky,  1998).  Abundant  <lpm 
bubble  walls,  sometimes  with  wavy,  unstable  forms,  are  consistent  with 
rupture  at  about  this  thickness  (Klug  and  Cashman,  1996).  Number-size 
distributions  for  vesicles  larger  than  a  few  tens  of  microns  have  power-law 
shapes  with  fractal  dimensions  of  — 1.9  to  —  3.4  attributed  to  coalescence 
(Gaonac’h  et  al.,  1996a,b;  Klug  et  al.,  2002).  Pumice  fractal  dimensions  are 
similar  to  those  at  Mount  Mazama  ( —  3.0  to  —  3.5;  Klug  et  al.,  2002),  and 
slightly  higher  than  those  at  Novarupta  (  —  3.9;  Adams  et  al.,  2006).  The 
degree  of  coalescence  increases  from  breadcrust-bomb  rinds,  to  bread- 
crust-bomb  interiors,  to  pumices,  as  witnessed  by  progressively  lower 
fractal  dimensions  (—1.9  to  —2.3,  —2.7  to  —2.8  and  —2.9  to  —3.4, 
respectively)  and  higher  proportions  of  intermediate-sized  vesicles  (0,  0 
to  9,  and  25  to  44%,  respectively).  This  is  attributed  to  the  progressively 
longer  time  interval  between  decompression  (and  bubble  nucleation) 
and  quench.  Moreover,  the  fraction  of  isolated  vesicles  is  in  general  lower 
in  pumices  than  in  breadcrust  bombs,  based  on  either  pycnometry  or 
image-analysis  data  (Table  1 ).  Advanced  coalescence  on  timescales  of  a 
few  to  a  few  tens  of  seconds  is  consistent  with  laboratory  experiments  on 
vesiculating  rhyolitic  melts  (Westrich  and  Eichelberger,  1994;  Navon  and 
Lyakhovsky,  1998;  Larsen  et  al.,  2004;  Burgisser  and  Gardner,  2005). 


Coalescence  continued  in  the  hot  interiors  of  >30  cm  pumices  for 
some  time  after  explosion,  forming  bubbles  up  to  a  cm  or  more  in  some 
metre-sized  clasts.  The  lack  of  systematic  radial  vesicularity  gradients 
and  surface  breadcrusting  show  that  this  must  have  taken  place  without 
expansion.  Coalescence  in  a  static  foam  is  driven  by  capillary  forces  and 
gravitational  drainage  of  melt  films  (Navon  and  Lyakhovsky,  1998).  It 
must  have  taken  place  on  a  timescale  shorter  than  that  required  for  the 
cores  of  large  pumices  to  cool  from  magmatic  temperature  (-850  °C) 
through  the  glass  transition  (-700  °C;  Dingwell,  1998).  The  time  for  the 
core  temperature  of  a  sphere  of  diameter  d  to  cool  conductively  by  20% 
in  contact  with  air  is  0.08d2/4/c,  where  k  is  thermal  diffusivity  (Carslaw 
and  Jaeger,  1959,  section  9.3.IV).  Thermal  diffusivity  of  pumice  is  a 
function  of  temperature,  with  a  value  of  3x  10“ 7  m2s-1  at  800  °C 
(Bagdassarov  et  al.,  1994);  internal  heat  transfer  by  convection  or 
radiation  can  be  neglected  (Hort  and  Gardner,  2000).  This  yields  -2  h 
for  a  30-cm  pumice,  which  is  consistent  with  the  authors’  observations 
(TF1D)  that  dm-sized  pumices  were  only  warm  to  the  touch  a  day  after 
emplacement.  The  value  for  a  1-m  pumice  is  -20  h,  so  that  static 
coalescence  may  have  continued  as  long  as  a  day  in  the  largest  pumices. 
These  estimates  ignore  premature  cooling  by  rain.  Any  air  circulation 
was  insufficient  to  oxidize  the  hot  clast  interiors  during  cooling  (Tait 
et  al.,  1998). 

8.6.  Eruptive  model 

The  results  are  now  placed  in  the  context  of  the  eruptive  model 
proposed  in  previous  papers  (Druitt  et  al.,  2002;  Melnik  and  Sparks, 
2002;  Clarke  et  al.,  2002;  Formenti  et  al.,  2003;  Diller  et  al.,  2006; 
Clarke  et  al.,  2007,  Burgisser  et  al.,  in  press)  and  illustrated  in  Fig.  15. 

Magma  ascended  by  viscous  flow  during  the  3-63  h  between 
explosions.  The  mean  ascent  velocity  during  this  time  was  >0.02  m  s-  \ 
corresponding  to  a  magmastatic  decompression  rate  of  0.001- 
0.005  MPa  s_1.  Many  large  vesicles  in  the  pumices  may  represent  bubbles 
nucleated  between  explosions,  the  number  densities  of  large  vesicles 
being  consistent  (via  Eq.  (1))  with  nucleation  at  such  decompression 
rates.  Other  large  vesicles  may  be  the  coalesced  remnants  of  bubbles 
nucleated  at  higher  number  densities  during  previous  explosions.  Inter¬ 
explosion  coalescence  generated  the  permeability  necessary  for  gas 
escape  from  the  ascending  magma,  as  reflected  in  the  low  water  contents 
of  groundmass  glasses  in  dense  pumices,  breadcrust  bombs  and  glassy 
lava  fragments  (Burgisser  et  al.,  in  press).  Possible  relict  gas-escape  paths 
are  preserved  in  dense  pumices  as  interconnected  networks  of  large, 
partially  collapsed  vesicles.  Alternating  bands  of  variable  vesicularity  in 
some  flow-banded  dense  pumices  and  finely  breadcrusted  bombs  may 
record  the  process  of  densification  induced  by  differential  outgassing, 


Table  2 

Estimation  of  timescales  pertaining  to  vesicle  growth  in  explosion  pyroclasts. 


a  Burgisser  et  al.  (in  press). 

b  Harford  et  al.  (2003).  The  temperature  of  the  breadcrust  bomb  rind  was  arbitrary  decreased  by  150  °C  to  take  into  account  the  quenching  by  the  atmosphere  under  the  glass 
transition  temperature. 

c  Method  ofToramaru  (2006),  using  Zhang  and  Behrens  (2000)  and  Behrens  et  al.  (2004),  with  Si02  contents  of  76  to  80  wt.%  for  all  samples  (Barclay  et  al.  1998),  cr=  0.025  N  m~ 1 
for  heterogeneous  nucleation  (Cluzel  et  al.,  2008)  and  the  given  melt  water  contents. 

d  From  Table  1 ,  Nv  values  of  both  BCP1  and  BCP43  used  for  the  breadcrust  bomb  core  and  rind,  values  of  AM036,  AM029  and  R2  for  the  poorly  vesicular  pumice  and  values  of  PV3 
for  the  highly  vesicular  pumice. 

e  Hess  and  Dingwell  (1996),  77  calculated  using  given  melt  water  contents. 

f  Calculated  following  Lensky  et  al.  (2004);  see  text  for  details. 
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Fig.  15.  Summary  of  vesiculation  processes  during  an  average  1997  Vulcanian  explosion. 


quenched  in  progress.  Immediately  prior  to  an  average  explosion,  the 
conduit  consisted  of  a  deep,  homogeneous,  gas-rich  zone  with  <55  vol.% 
voidage  overlain  by  an  impermeable,  degassed,  microlite-rich  and 
rheologically  stiffened  plug,  behind  which  gas  pressure  accumulated  in 
a  transition  region  of  partly  degassed  magma  (Diller  et  al.,  2006;  Burgisser 
et  al.,  in  press;  Fig.  15a). 

Plug  rupture  occurred  once  overpressure  exceeded  a  critical 
threshold,  probably  of  a  few  MPa  or  more  (Druitt  et  al.,  2002;  Spieler 
et  al.,  2004b;  Burgisser  et  al.,  in  press).  Decompression  caused  upward 
expansion,  acceleration,  fragmentation  and  eruption  of  the  conduit 
contents  (Melnik  and  Sparks,  2002;  Mason  et  al.,  2006).  Coarsely 
breadcrusted  bombs,  finely  breadcrusted  bombs,  dense  pumices  and 
light  pumices  form  a  continuous  assemblage,  the  textural  features  of 
which  can  be  understood  by  differences  in  the  relative  timing  of 
initial  decompression,  fragmentation,  and  surface  and  interior 
quench.  Largely  degassed,  non-vesicular,  microlite-rich  magma  of 
the  upper  part  of  the  plug  underwent  decompressions  of  a  few  MPa 
and  entered  the  atmosphere  immediately  after  fragmentation.  Cool¬ 
ing  during  flight  limited  bubble  growth  in  surface  layers,  whereas  the 
hotter  interiors  frothed  up,  cracking  the  outer  rind  and  breaking 
amphibole  phenocrysts.  These  clasts,  characterized  by  a  large  vesi- 
cularity  difference  between  core  and  rind,  form  the  coarsely  bread¬ 
crusted  bombs  (Fig.  15b).  Dense  magma  fragments  from  deeper  in  the 
plug  followed  a  similar  history,  but  the  time  interval  between  decom¬ 
pression  and  discharge  was  slightly  longer,  allowing  more  bubble 
growth  in  the  surface  rind  prior  to  quench.  The  vesicularity  contrast 
between  surfaces  and  cores  of  the  finely  breadcrusted  bombs  is 
therefore  smaller.  Some  coarsely  breadcrusted  bombs  that  broke 
either  during  ballistic  impact  or  transport  in  pyroclastic  flows  de¬ 
veloped  finely  breadcrusted  surfaces  on  breakage  surfaces.  In  all 
bombs,  core  inflation  and  surface  breadcrusting  was  complete  long 
before  pyroclastic  flows  came  to  rest. 

The  time  interval  between  decompression  and  discharge  of 
the  deep,  gas-rich  magma  was  sufficiently  long  that  syn-explosive 
inflation  was  complete  prior  to  surface  quench  (Fig.  15c).  This 
accounts  for  the  lack  of  surface  breadcrusting,  the  apparent  lack  of 
dependence  of  vesicle  textures  on  clast  size  (up  to  ~30  cm),  and  the 


textural  similarities  between  flow  and  fallout  pumices,  despite 
different  thermal  histories  after  leaving  the  vent.  It  also  allowed 
coalescence  of  newly  formed  bubbles  to  proceed  further  than  in  the 
breadcrust  bombs.  As  pumices  inflated,  many  amphibole  phenocrysts 
were  broken  in  a  manner  suggestive  of  isotropic  decompression  in 
parcels  of  fragmented  magma.  However  uniaxial,  flow-foliation- 
parallel  boudinage  of  amphiboles  implies  that  syn-explosive  vesicu¬ 
lation  of  many  pumices,  possibly  those  from  deeper  in  the  conduit, 
was  essentially  complete  even  prior  to  fragmentation. 

Syn-explosive  bubble  nucleation  at  all  levels  of  the  conduit  took 
place  heterogeneously  on  titanomagnetite.  Growth  took  place  under 
disequilibrium  conditions  largely  by  decompression-driven  expan¬ 
sion;  diffusive  outgassing  of  melt  was  limited  by  the  short  timescale  of 
the  explosions,  accounting  for  the  high  water  contents  preserved  in 
groundmass  glasses. 

Decompression  rates  estimated  from  heterogeneously  nucleated 
vesicle  number  densities  (at  least  0.3-6.5  MPa  s-1)  are  an  order  of 
magnitude  lower  than  if  homogeneous  nucleation  were  assumed 
(e.g.  Toramaru,  2006,  Hamada  et  al.,  in  press).  However  they  may  be 
under-estimated  if  there  was  insufficient  time  for  bubble  numbers  to 
reach  stationary  values.  Plug  rupture  during  Vulcanian  eruptions 
generates  a  downward-propagating  rarefaction  wave,  followed  by  a 
slower  fragmentation  front  (Melnik  and  Sparks,  2002),  the  ascent 
speed  of  the  bubbly  magma  being  a  great  deal  slower  than  either. 
A  zone  of  steep  pressure  gradient  develops  beneath  the  fragmentation 
front  due  to  the  discontinuous  pressure  change  at  the  fragmentation 
surface  and  wall  friction  (Koyaguchi  et  al.,  2008).  Individual  parcels  of 
magma  in  an  average  Montserrat  explosion  simulated  by  Melnik  and 
Sparks  (2002;  their  Fig.  4)  experience  decompression  rates  of  ~1  to 
>20  MPas-1  as  they  approach,  then  pass  through,  the  descending 
fragmentation  front.  These  values  are  of  the  same  order  of  magnitude 
as  those  estimated  from  vesicle  number  densities,  showing  that  there 
is  no  incompatibility  between  observations  and  models.  They  are 
consistent  with  the  idea  that  syn-explosive  bubble  nucleation  during 
the  1997  explosions  was  initiated  during  rapid  decompression  of 
conduit  magma  in  a  zone  of  steep  pressure  gradient  beneath  the 
descending  fragmentation  front. 


230 


T.  Giachetti  et  al.  /  Journal  of  Volcanology  and  Geothermal  Research  193  (2010)  215-231 


Acknowledgments 

The  project  was  partly  funded  by  the  French  Agence  Nationale  de 
la  Recherche  (ANR-EXPLANT,  Contract  No  ANR-05-CATT-003  to  C. 
Martel).  We  benefited  from  discussions  with  Stephane  Poussineau, 
Nicolas  Cluzel,  Karim  Kelfoun,  Thomas  Shea,  Lucia  Gurioli,  Didier 
Laporte,  Jean-Louis  Bourdier  and  Jim  Gardner.  Marc  Leblanc  gave  us 
access  to  the  He  pycnometer,  and  Jean-Marc  Henot  helped  us  use  the 
SEM.  We  thank  Amanda  Clarke  and  Oliver  Spieler  for  their  thoughtful 
and  constructive  reviews  of  the  manuscript. 

Appendix  A.  Size-distribution  measurement  of  vesicles  and  crystals 

Size  distributions  were  measured  using  images  of  six  overlapping 
magnifications:  (1)2  images  of  a  thin  section  and/or  cut  hand  specimen, 
(2)  2  SEM  images  at  x27,  (3)  2  SEM  images  at  x  100,  (4)  4  SEM  images 
at  x200,  (5)  4  SEM  images  at  x500,  and  (6)  6  SEM  images  at  xl500. 
Vesicles  or  crystals  (referred  to  hereon  as  ‘objects’)  of  >300  pm  were 
counted  at  magnification  1,  and  those  of  <300  pm  were  measured  at 
magnifications  2-6.  Size  distributions  were  constructed  in  four  steps 
as  follows: 

Step  1.  Size  distribution  for  each  image.  Each  image  resulted  in 
three  binary  images  (isolated  vesicles,  connected  vesicles,  and 
crystals),  which  were  analysed  with  the  ImageJ  image-analysis 
software.  Using  the  surface  area  A  and  the  scale  of  the  image,  an 
equivalent  diameter  D  was  calculated  for  each  object  (assuming 
circular  shape),  and  surface  area  and  area  number  distributions 
were  calculated  as  functions  of  D.  The  probability  of  a  large  object 
being  intersected  by  the  edge  of  the  image  is  larger  than  that  of  a 
small  object,  and  we  applied  the  corrections  of  Miles  (1974)  and 
Lantuejoul  (1980)  to  obtain  the  true  number  of  objects  (Na). 

Step  2.  Size  distribution  for  each  magnification.  The  object  size 
distribution  for  a  given  magnification  was  obtained  by  averaging 
those  of  the  separate  images  for  that  magnification.  This  resulted 
in  averaging  any  heterogeneity  in  object  size  or  shape  distribution 
present  on  the  microscopic  scale,  assuming  that  possible  shape 
variations  were  randomly  oriented.  Tests  on  multiple  areas  of 
single  pyroclasts  showed  that  this  assumption  was  valid  to  a  first 
approximation. 

Step  3.  Size  distribution  for  the  whole  sample.  This  third  stage 
generated  an  object-size  distribution  for  the  whole  sample  using  the 
individual  distributions  at  each  magnification.  Some  studies  perform 
this  step  using  Na  (Shea  et  al.  2010;  Klug  and  Cashman,  1994;  Klug 
et  al.,  2002),  whereas  we  used  object  surface  areas.  We  first  defined 
cut-off  limits  between  magnifications,  being  chosen  so  that  each 
object  had  a  diameter  of  at  least  1 5  pixels  and  was  no  larger  than  1 0% 
of  the  image  size  (in  order  to  be  statistically  represented).  The  limits 
correspond  to  log  D  values  of  2.5,  2.0, 1.6, 1.2  and  0.8,  from  lowest  to 
highest  magnifications,  as  shown  on  Fig.  10.  We  then  normalised  the 
data  at  each  magnification  to  the  surface  area  of  the  lowest 
magnification,  correcting  at  each  magnification  for  objects  larger 
than  the  cut-off  limits  previously  defined.  This  correction  was  applied 
starting  at  the  lowest  magnification  and  working  progressively 
upwards.  Surface  and  number  distributions  were  then  obtained  for 
the  whole  sample.  In  total  930-4750  crystals  and  3390-9620  vesicles 
were  measured  for  each  sample.  Full  explanation  and  worked 
examples  are  available  from  the  first  author  on  request. 

Step  4.  2D-3D  conversions.  Two-dimensional  data  (Na  and  surface 
distribution)  were  then  converted  to  3D  ( Nv  and  volume  distribu¬ 
tion).  Since  we  chose  samples  with  spherical  and  homogeneously 
distributed  vesicles  (at  least  to  the  eye),  we  assumed  that  the  area 


distribution  obtained  directly  equals  the  volume  distribution.  Our 
method  has  the  advantage  of  directly  providing  the  total  sample 
vesicularity  by  integrating  the  vesicle-volume  distribution.  Values 
thus  obtained  were  close  to  those  from  He-pycnometry  (Table  1), 
convincing  us  of  the  accuracy  of  our  method.  Two  methods  were 
used  to  convert  Na  to  Nv.  Cheng  and  Lemlich  (1983)  proposed  a 
conversion:  Nv  =  Na/Dm,  where  Dm  is  the  average  equivalent 
diameter  of  each  size  class.  This  method  is  easy  to  apply  and 
involves  no  assumptions  of  object  shape.  However,  it  does  not  take 
into  account  the  cut  effect,  which  is  corrected  by  the  method  of 
Sahagian  and  Proussevitch  (1998),  assuming  spherical  shapes.  Our 
results  showed  that  object  number  densities  obtained  from  the  two 
methods  were  the  same  to  within  15%,  and  mostly  within  5% 
(Table  1).  Number  density  distributions  were  also  very  similar. 
However,  since  the  method  of  Sahagian  and  Proussevitch  (1998) 
sometimes  generates  negative  values  of  Nv  for  some  size  classes, 
values  used  in  the  text  are  those  obtained  using  the  method  of  Cheng 
and  Lemlich  (1983). 

Appendix  B.  Supplementary  data 

Supplementary  data  associated  with  this  article  can  be  found,  in 

the  online  version,  at  doi:  10.1 01 6/j.jvolgeores.201 0.04.001. 
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